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In this paper, a generalized three-phase multilevel power inverter (MLI) structure is 
proposed under asymmetric configurations. The operating mode and the switching 


combinations are briefly exposed according to the parity of the number of direct cur- 
rent (DC) voltage sources in use. Subsequently, the proposed topology is evaluated in 
terms of commonly used factors and then benchmarked against some of the state-of- 
the-art cascaded MLIs featuring multiple DC voltage sources (MDCS-CMLIs) while 
putting emphasis on the reduction of power switching devices. Moreover, a new near- 
est level control (NLC)-based modulation technique is designed for the purpose of 
better comply with some quality grid codes, namely the European EN 50160 and the 
International IEC 61000-2-12. The identification of the optimal control thresholds is 
realized by a constrained optimization algorithm (e.g., particle swarm optimization 
(PSO)) which is implemented in python script and validated through SIMULINK fast 
fourier transform (FFT) analysis tool. Lastly, the harmonic performance of the pro- 
posed technique is compared side-by-side with that of the conventional NLC scheme 
and exhibits significant reduction in harmonic distortion. 
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1. INTRODUCTION 

In medium-voltage (MV) and medium-to high-power level conversion field, multilevel inverters (MLIs) 
are considered to be the flagship solution for grid-related applications including renewable energy sources in- 
tegration, high voltage DC transmission/flexible AC transmission system (HVDC/FACTS), shunt active power 
filter (APF), and static synchronous compensator (STATCOM). Moreover, they are renowned as one of the most 
effective way to synthesize high resolution AC power endowed with attractive features such as low harmonic 
distortion, high power handling, and best power conversion efficiency to name a few [1]. Despite all of the 
aforementioned benefits, MLI topologies require a relatively large number of power semiconductor and passive 
components. To address this issue, a lot of research are led toward the reduction of the power components 
count so as to enhance both the overall power efficiency and the system reliability [2]-[4]. 

In practice, cascaded MLI structures with multiple DC sources (MDCS-CMLIs) are regarded to be 
the most authoritative MLIs since they have been employed in some renewable energy systems for more than a 
decade. They can be set either in symmetric or in asymmetric configuration [1], [5]. Asymmetric nature of DC 
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sources would bring more levels to the output waveforms which positively impact their harmonic performance. 
For this reason, only asymmetric structures are considered throughout this paper. 

Low switching frequency modulations (LSFM) found their way in the field of MLI structures with the 
main objective of alleviating the switching losses incurred by the relative high number of power switches [6], 
[7]. Nowadays, there is only one well-established LSFM scheme geared toward grid-tied applications known 
as selective harmonic mitigation (SHM) which in turn is akin to the selective harmonic elimination (SHE) 
modulation technique [8]-[11]. However, SHE/SHM has some limitations that could be overcome by NLC 
[5], [12], especially those related to dynamic performance under real-time operation and computational effort 
required while using lookup tables or curve fitting polynomials or any initial guessing and iterative routines 
[13], [14]. 

The next section proposes a generalized three-phase MLI topology under different parity of k DC 
voltage sources. To evaluate performance, the inverter structure is analyzed in terms of some commonly used 
factors and compared to other state-of-the-art topologies. Afterwards, in section 3, a new NLC-based mod- 
ulation is developed for the purpose of better meeting the grid quality requirements of both EN 50160 and 
IEC 61000-2-12 codes. In order to achieve that, some mathematical formulations based on fourier series ex- 
pansion are shortly exposed. Next, the formulation of the optimization problem is given as a system of nonlinear 
inequalities considering each individual harmonic limit L, and the total harmonic distortion (THD) of a given 
line-to-line output voltage ve(wt). Then, the particle swarm optimization (PSO) algorithm [15] is called to 
determine the appropriate switching angles 6; depending on the voltage ratio r and the thresholds A;. Finally, 
the harmonic performance of the proposed control technique is benchmarked side-by-side against that of the 
conventional NLC scheme. Results and discussion are carried out while relevant conclusions are highlighted at 
the end of this paper. 


2. ASSESSMENT OF THE PROPOSED MLI TOPOLOGY 
2.1. Basic unit cell MDCS-CMLI structure 

The proposed three-phase MLI topology is depicted by Figure 1 where the basic unit cell is shaded in 
light color. As it can be seen, each single-phase cell is made up from two voltage sources, four unidirectional 
(511, 81,2, 53,1, 53,2) and two bidirectional power switches ($31. S39). Several unit cells can be strung to- 
gether in series so as to synthesize high resolution output voltages. A comparative study with detailed operating 
modes for both 5-level symmetric and 7-level asymmetric configurations has already done in [16], so it will not 
be discussed further except for some reminiscent results required when designing the generalized structure. 


CELL #1 CELL #2 see CELL #m 


THREE-PHASE 
LOAD 


Figure 1. Proposed generalized three-phase multilevel power inverter topology. Basic unit cell is shaded in 
which bidirectional switches are illustrated by two back-to-back N-channel power MOSFETs 


The focus of the present work was mainly driven toward grid-connected MLIs for renewable energy 
integration. Asymmetric nature of the discrete DC sources would be of great benefit in enhancing power quality 
required by the grid code. In each cell, a binary voltage ratio was used (i.e., r = V2/V, = 2) and for the sake 
of efficiency, bidirectional switches were made of two back-to-back N-channel power MOSFETs. 

Two aspects of importance should be at the forefront when designing MLI topologies: total standing 
voltage (TSV) and total active switches (TAS). The former is defined as the sum of the highest voltage stress 
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across each power switch while the latter is the total number of conducting MOSFETs over one time period: 


4(V, + Vo) + 2V2 = 16V; (1) 
32 MOSFETs. 


TSV 
TAS 


2.2. Extended MDCS-CMLI structure 


The basic unit cell within the generalized structure of Figure | can be extended in two different ways 
according to the parity of the k number of DC voltage sources. For each phase, every unit cell is connected 
together in series so as to generate high-resolution output waveforms by virtue of arithmetic combinations of 
the DC power sources where their p.u. magnitudes Vy feature septenary configurations as follows: 


7k-D/2, k=2m-1 


Ym € N*, Vk/ Vi = 
ʻ k/Vi = 1 2(76/2-1), & = 2m, 


(2) 


where m represents the number of basic unit cells being connected. 
In order to build up all the available n levels and guarantee every time that any of the 2m DC voltage 
sources will never be short-circuited, then the following general rule should be applied: 


2m 3 


S= {0/1} D Sg 1, (3) 


j=1 i=1 


Si j correspond to the enhancement-mode NMOS, so ‘0’ and ‘1’ denote, respectively, the OFF- and ON-state. 
In addition, the available switch combinations for the generalized single-phase MLI are tabulated in Table 1. 


Table 1. Switching combinations of the generalized single-phase MLI under asymmetric operation 


PER-UNIT OUTPUT VOLTAGE LEVELS (Vk /V1) 


1-7" 7m 4 
ve) —24 --- -10 --- -3 -2 -1 O 4l 42 43 --- 410 ++. 424 -:- 
2 2 

Sit ON > ON --- ON --- ON OFF ON ON OFF OFF OFF --- OFF --- OFF --- OFF 
S12 OFF --- OFF --- OFF --- OFF OFF OFF ON ON OFF ON --- ON --- ON -> ON 
S13 ON > ON- ON --- ON ON ON ON ON ON ON -:- OFF --- OFF + OFF 
Si OFF -:- OFF --- OFF --- ON ON ON ON ON ON ON -:- ON -:- ON -> ON 
Si2m-1 ON > ON- ON --- ON ON ON ON ON ON ON- ON- ON -> OFF 
S1,2m OFF > ON >- ON --- ON ON ON ON ON ON ON ON- ON -> ON 
Sž OFF --- OFF --- OFF --- OFF ON OFF OFF ON OFF OFF --- OFF --- OFF --- OFF 
Gis OFF --- OFF --- OFF --- OFF OFF ON OFF OFF ON OFF --- OFF --- OFF -:: OFF 
cn OFF --- OFF --- OFF --- OFF OFF OFF OFF OFF OFF OFF --- ON --- OFF --- OFF 
Sii OFF --- OFF --- ON --- OFF OFF OFF OFF OFF OFF OFF -:- OFF --- OFF --- OFF 
S$om—1 OFF ++: OFF + OFF --- OFF OFF OFF OFF OFF OFF OFF --- OFF -+ OFF --- OFF 
Cn OFF --- OFF --- OFF --- OFF OFF OFF OFF OFF OFF OFF -:- OFF --- OFF --- OFF 
S31 OFF --- OFF --- OFF --- OFF OFF OFF OFF OFF ON ON >- ON --- ON -> ON 
S32 ON > ON --- ON --- ON ON OFF OFF OFF OFF OFF --- OFF --- OFF --- OFF 
53.3 OFF --- OFF --- OFF --- OFF OFF OFF OFF OFF OFF OFF --- OFF --- ON -= ON 
S34 ON > ON --- OFF --- OFF OFF OFF OFF OFF OFF OFF --- OFF --- OFF --- OFF 
S3,2m-1 OFF -:- OFF --- OFF --- OFF OFF OFF OFF OFF OFF OFF --- OFF --: OFF ON 
S'3,2m ON .-- OFF --- OFF --- OFF OFF OFF OFF OFF OFF OFF --- OFF --- OFF OFF 


Depending on the parity of k number of DC voltage sources, the equations of the number of switches 
(Nsw), gate drivers (Nga) and total standing voltage (TSV) for the proposed n-level power inverter are given 
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as shown in: 
Ym € N* | k = 2m, Ym € N* | k= 2m- 1, 
n = 72 n = 3 (7-9/2) 
Nsw = 8log;(n) Nsw = 8log;(n/3) +4 4) 
Noa = 6log;(n) Nga = 6logz(n/3) +4 
TSV = 5$; 16V;/3 TSV = 5# 16V;/3 + 4V 
= 8(n — 1)V1/3. = (20n — 24) V; /9. 


It is worth noting that the use of the last two bidirectional power switches (i.e., S3 sm-1 and S3 om> 
m > 2) is not justified when k is odd; only unidirectional are needed. Furthermore, and in order to estimate 
the overall cost-per-level involved in creating an n-level inverter structure, the cost function $(n, 3) is taken 
into consideration with the assumption that all the switches in use share the same power rating [17], [18]. 
Depending on its values (greater or less than unity), the weight coefficient 8 is used to emphasize either the 
TSV pu or the total number of required electric components: 


$(n, B) = (New + Nga + [Ne + Na] + BTSVpu) n! k, (5) 


where the per-unit TSV is defined as: TSV pu = TSV/ oF V;. However, topologies with power diodes N4 
and/or electrolytic capacitors N are not considered here since they are generally lossy and/or bulky. 


2.3. Comparative study with recent MDCS-CMLIs 

A comparative study of some cutting-edge MDCS-CMLIs is performed in order to show the merits 
of the proposed structure [19]-[22]. Typically, the n is limited to 50 levels, which would be far enough for 
any grid-tied applications targeting either single- or three-phase inverters. As can be seen from Figure 2, the 
proposed topology requires the least number of k(n) DC sources with moderate Nsw(n), Nga(n) and p.u. 
TSV (n), resulting in the lowest cost function $(n, 8) for both cases (i.e, 8 < 1 and 6 > 1) regardless the 
parity of k. 


—e— Proposed (k even) —#— Proposed (k odd) [19] [20] —— [21] —— [22] 


Figure 2. Comparison of the proposed asymmetric topology to some state-of-the-art MDCS-CMLIs 


3. PROPOSED CONTROL STRATEGIES FOR GRID-TIED MLI TOPOLOGIES 
3.1. Mathematical formulation of the proposed MLI output voltage waveform 

The output voltage waveform of the proposed inverter driven by an NLC modulation is similar to a 
staircase sine wave function with equal riser (r = 2), and because of its odd quarter-wave symmetry property, 
the general form of the Fourier coefficient b, is simplified as shown in with ap = b2} = 0 for all harmonics A: 


4 
m==( 


T 


02 03 04 
i) sinapa + f rsinhodd + f (r + 1) sin hé dé 
Oy 02 03 


n/2 log (n) (6) 
toot f (r +1) 5 (T"= +1) sinh a8) 


On-1 


5 m=2 
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The fourier series of the phase voltage vg(wt) can then be written as a function of the voltage ratio r 
and the switching angles 6;: 


vg(wt) = », bn sin hut = > [cos h01 + (r — 1) cos hO2 + cos hO3 + (5 — 2r) cos h04 + cos hOs 
+ (r — 1) coshOg +--+ + cos hôin—1)/2] Sin hwt (h : odd harmonics) 
where 0; is defined by (8), taking into account the number of available levels n and the thresholds \;. The 
introduction of the latter parameters allows to distinguish two operating modes regarding the symmetry or 


asymmetry of the thresholds in which the magnitudes of A; are, respectively, equal or not necessarily equal. 
Thus, conventional NLC modulation can be re-established simply by setting all symmetrical thresholds to 1. 


2\;(i — a 2 


Vi € {1,2,...,(n—1)/2} JA € RỌ YA; O< Oh < Bay. 65 = arcsin ( i = (8) 
"E 


The generalized form of the modulation index m,(r, A) is given is being as: 


a Îsi _ 4,17, A) 
mas) max(igi)  2(n +2r-— 5)Vı 0) 


where 64 ,1(r, A) and max(ôg 1) represent, respectively, the amplitude of the fundamental component of vg (wt) 
and its maximum by which all 6; are zeroed out. 


3.2. Constrained optimization problem for grid-tied applications 

Power quality grid codes usually specify an upper limit Lp for each specific harmonic h and indicate 
the maximum accepted distortion which, conventionally, considers all harmonics up to the 40th. For instance, 
the signal quality requirements of both the European and the International grid codes, namely EN 50160 and 
IEC 61000-2-12, are tabulated in Table 2, where the limited voltage THD is set to 8%. It is worth noting that 
the maximum harmonic orders required by EN 50160 and IEC 61000-2-12 are 25th and 50th, respectively [23], 
[24]. 


Table 2. Individual harmonic limits according to codes EN 50160 and IEC 61000-2-12 in MV networks 


EN 50160 (h <25) IEC 61000-2-12 (h <50) 

ODD HARMONICS ODD HARMONICS EVEN HARMONICS ODD HARMONICS ODD HARMONICS EVEN HARMONICS 
(Non-multiple of 3) (Multiple of 3) (Non-multiple of 3) (Multiple of 3) 
h Ln (%) h Ln (%) h Ln(%) h Ln (%) h Ln (%) h Ln (%) 
5 6 3 5 2 25 6 3 5 2 2 
7 5 9 1.5 4 17 5 9 1.5 4 1 
11 3.5 15 0.5 >6 0.5 11 3:3 15 0.4 6 0.5 
13 3 21 0.5 13 3 21 0.3 8 0.5 
17 2 17 2 >27 0.2 10 0.5 
>19 1.5 > 19 38.59/h — 0.27 >12 2.5/h +0.25 


To seek the appropriate voltage ratio r and all the thresholds A; (i.e., the switching angles @;) that 
would have compliance with the aforementioned quality grid regulations, a constrained optimization problem 
whose cost function CF (r, A) is formulated is being as: 


CF(r, A) = CF(r, A1, 2,3, -- <, Am—1)/2) = min THD.o [ve(wt)] 


2v 3V) 
E E E ) 
T 
AV/3V, 
s.t. (1 =a Ô)Ln 2 impaled) [cos ho, + (r =e 1) COS hos + COS h03 + RIA + Cos hO(n—1)/2] 
37 ^2 
ô p (r, A) 
> THD t) = > > 
A 7 . ne ) h=5,7 ôg (r, A) 


(10) 
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where 6 is a safety margin; ve and vg p are, respectively, the line-to-line output voltage and its specific non- 
triplen odd harmonics under three-phase balanced condition. They are interrelated and can be written as shown 


in: 


ve(wt) = V3 5 bn cos(hwt + phr) = 5 ve nlr, A) (h : non-triplen odd harmonics) (11) 
h=1 h=1 


3.3. Assessment and implementation of the proposed modulation technique 

MLIs with higher number of levels are more likely to meet grid quality requirements. However, they 
significantly suffer from low power efficiency, bulky footprint, high cost and complex design and control. 
For instance, the proposed inverter structure with 21 levels may be employed, but instead it would be more 
advisable to use that with 7 levels tuned by a small passive filter since the number of power switches N.w(n), 
DC sources k(n) and TSV (n) are reduced by factors of 1.75, 2, and 2.75, respectively. Figure 3 depicts the 
integration of the proposed 7-level power inverter into the medium-voltage (MV) power distribution system. 


POWER INVERTER 


GRID 
| LINE FREQ. PCC (20/22 kV) 
seat TRANSFO. Lg 


Ox ge iaa 
FNAS 
al Y A Sw 


Figure 3. Integration of the proposed 7-level power inverter into medium-voltage power distribution networks 


LOAD (Z,) 
š 


Henceforth, all the above generalized equations are called to deal with only one symmetrical threshold 
under A-NLCg modulation or with three asymmetrical thresholds under A-NLC 4. The domain of each of these 
thresholds depend upon the definition given by (8), where in the first case A € [0,1.2] and in the second 
case Ay € [0,6], Ag € [0,2] and As € [0,1.2] such that Ay < 3A < 5A3. The major hindrance of solving 
the constrained optimization problem arises from the challenge of exploring a wide search space in order to 
locate feasible solutions {\*,7r*} that give best harmonic profiles of the output voltage. So, with the intention 
to enhance the solutions quality and speed up the processing time of the used evolutionary algorithm, only 
potential feasible solutions with THD. < 8% are kept and then plotted in Figures 4 (a) and (b) for both à- 
NLCg and A-NLC 4 modulation strategies, respectively. Moreover, it can be seen in Figure 4 (c) that MLIs 
driven by conventional NLC (A = 1) are out of the allowed range and those with binary configuration (r = 2) 
exhibit the lowest THD. The last result has already been demonstrated in [25]. 


4.0 
3.5 
THD > 8% = 
3.0 Bs 
= S 
Sas 2 
i E 
m 4 
2.04 E 
1.54 
1.0 
0.4 
5.0 5.5 6.0 6.5 7.0 7.5 8.0 
THD% (%) 
(a) (b) (c) 


Figure 4. Potential feasible solutions to the constrained optimization problem under: (a) A-NLCg and (b) 
A-NLC 4 (with r = 2) modulations. (c) Best harmonic performance with respect to voltage ratio r 
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In consequence, the binary configuration is chosen which additionally narrow down the search space. 
The best possible solutions {A} to the non-linear constrained optimization problem of (11) are determined 
using an implementation of PSO in Python. The validation is done by the algorithm of Figure 5 under 
SIMULINK. Figure 6 plots the modulating signals, phase- and line-output voltages under conventional, sym- 
metrical and asymmetrical NLC modulations, whereas Figure 7 indicates, side-by-side, the harmonic contents 
of the respective line-output voltages along with the limits imposed by European and International Grid Codes. 
Noticeable enhancements in terms of harmonic distortion are directly obtained by the proposed A-NLC s and 
A-NLC 4 modulations, particularly the THD reductions are of about 34% and 43% with respect to the conven- 
tional NLC when targeting full EN 50160 compliance with 10% of safety margin, and about 30% and 41% 
when targeting IEC 61000-2-12 Grid Code. It is worth to mention that all the lower order harmonics are can- 
celed out whereas few of those with higher order can be mitigated by introducing small low-pass filters. These 
results are summarized within Table 3 in which best feasible solutions are recorded. 


IN : Mod. signal SIG(¢) = (1 + r)Vı sin(wt + y) 
OUT : Gating signals for power switches S; j case equals to 1 
| Sij 4+ switch state #+1 (cf. Table 1); break 


if SIG(t) < AV and SIG(t) > —à Vi then case equals to 2 
a a | Si; — switch state #+2 (cf. Table 1); break 


else 
if SIG(t) > 2A2Vi or SIG(t) < —2X2V, then 
A © 3 return 


otherwise 


else 
if ROUND(SIG(t)/(AVi)) > +3 then 
[Ea ee | Sij < switch state #43 (cf. Table 1) 
switch ROUND(SIG(t)/(AV;)) do if ROUND(SIG(t)/(AV,)) < —3 then 
| Si j <- switch state #—3 (cf. Table 1) 


Figure 5. NLC-based algorithm for the implementation of the proposed 7-level modulation schemes 


NLC (ma = 0.80) THD19 = 8.8% A-NLCsg (0.95) THDao = 6.2% A-NLC4 (0.93) THDao = 5.1% 
F 3 +550 F 7 +550 F 5 


—550 E t i + —550 & = 
30 40 0 10 30 40 0 10 30 40 


0 10 


20 20 
TIME (ms) TIME (ms) 


20 
TIME (ms) 


— Modulating signal —— Line-to-Neutral output voltage vg(wt) ——— Line-to-Line output voltage ue (wt) 


Figure 6. Modulating signals, phase-and line-output voltages using NLC, A-NLCg and A-NLC, modulations 
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Figure 7. Harmonic contents of the line-output voltages using NLC, A-NLCs and A-NLC 4 modulations 
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Table 3. Salient results of the proposed modulation schemes under some best feasible solutions 


MODULATION THRESHOLDS THD GRID CODE COMPLIANCE NOTES 
SCHEME Ni XS A3 (%) EN 50160 IEC 61000-2-12 (Remaining harmonics) 
Conv. NLC 1.00 1.00 1.00 881 No No 13th-19th, 25th, 29th, 35th—41st 
\-NLCg 0.55 0.55 0.55 5.83 Yest N/A t Safety margin ô = 10% 

0.52 0.52 0.52 6.17 Yes No 31st, 35th and 37th 
\-NLC 4 0.61 0.56 0.68 5.01 Yest N/A t Safety margin ô = 10% 

0.52 0.55 0.68 5.15 Yes No 35th and 37th 


4. CONCLUSION 


In this article, a generalized three-phase cascaded multilevel inverter topology with reduced power 
switching devices was presented. The asymmetric nature of its odd/even DC power sources featuring a septe- 
nary relationship in voltage magnitude allows to synthesize high-resolution output waveforms. Subsequently, 
comparison to some of the cutting-edge MDCS-CMLIs were carried out. It has shown that the proposed struc- 
ture accomplishes the lowest cost-per-level within all levels regardless the value of the weight coefficient 2. 
Moreover, a new NLC-based modulation strategy was suggested with a twofold aim: firstly, to comply with 
two European and International power quality grid codes, respectively, EN 50160 and IEC 61000-2-12; and 
secondly, to provide an alternative pathway to overcome the intrinsic drawbacks inherited from SHE/SHM 
techniques. The obtained results within SIMULINK FFT analysis tool clearly demonstrated the outperfor- 
mance of the proposed modulation with (a) symmetrical thresholds A; over the conventional NLC, in which 
the optimal feasible solutions Až to the constrained optimization problem was determined by a bio-inspired 
algorithm (e.g., PSO) described in python programming language. 
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